Spatial Entanglement From Off-Diagonal Long Range Order in a BEC 
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We investigate spatial entanglement - particle-number entanglement between regions of space 
- in an ideal Bosonic gas. We quantify the amount spatial entanglement around the transition 
temperature for condensation (Tc) by probing the gas with two localised two-level systems. We 
show that spatial entanglement in the gas is directly related to filling of the ground state energy 
level and therefore to the off diagonal long range order (ODLRO) of the system and the onset of 
condensation. 
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I. INTRODUCTION 

Entanglement has been the subject of intense investi- 
gation since it was realised that its presence in quantum 
physics allows for new fast algorithms in computation 
For this, coherent control over a large number of 
qubits is needed, so it makes sense to study systems that 
already possess an amount of natural entanglement, for 
instance spin systems Q and between spatial regions in 
a Bose gas 0, [j] , then one can be sure that under cer- 
tain conditions entanglement exists even if the system is 
open. However, the entanglement found in such systems 
is more than just a computational resource, it is related 
to phase transitions in spin chains @, M, Q and to super- 
conductivity and superfluidity [Io| . 

Normally, many-body entanglement in such systems is 
detected, sometimes using macroscopic variables, such as 
temperature [TlJ] , magnetic susceptibility [l2| and heat 
capacity [l3| . that witness the entanglement. One can 
then extract the entanglement for other purposes [l4[ 
using probe systems. 

In this paper, unlike these previous approaches, we 
use two localised probe systems to measure the entangle- 
ment between spatial regions in a Bose gas. We refer to 
this type of entanglement as spatial entanglement, which 
manifests itself as non-local particle number correlations 
between different regions of space. Therefore, by locally 
probing these regions with effective qubits, we skip the 
need to find an entanglement witness. We can investi- 
gate entanglement without introducing probes [H , but in 
order to make our procedure fully operational we want to 
show that the entanglement is not an artefact of mathe- 
matics. On qubit-type systems standard experiments for 
the existence of entanglement can be run thereby sim- 
plifying the discussion of entanglement in a continuous 
variable (CV) system. 

In this paper we will show explicitly that spatial en- 
tanglement only exists between two distant localised vol- 
umes in a Bosonic gas below the critical temperature for 
a BEC (Tc). Thus the existence of spatial entanglement 
is directly related to the off-diagonal long range order 
(ODLRO) of the system. In fact it is the order in the 
Bose gas below Tc that allows spatial entanglement be- 
tween two distant regions to form. 

In the remaining parts of the introduction we will 



discuss the concept of entanglement between regions of 
space in II Al and discuss the notion of ODLRO in II Bl 
The rest of the paper will be dedicated to calculating the 
spatial entanglement contained in a Bose gas. Firstly we 
discuss the interaction of the probes and the gas in III Al 
followed by calculation of the amount of entanglement 
between the probes in sections III Bl and III CI Finally the 
results will be discussed in section III Dl and in the con- 
clusions in IIIII 



A. Entanglement between spatial modes 

Let us take a closer look at spatial entanglement it- 
self. When investigating entanglement between spatial 
regions in a CV system, one must work in second quan- 
tisation, where one speaks about the occupation num- 
ber of modes, because to correctly define entanglement 
the Hilbert space of a given system must have a tensor 
product structure. In the first quantisation the Hilbert 
space of two spatial regions is a direct sum of each region 
and therefore does not have the required tensor product 
structure. However, in this paper we use second quanti- 
sation where two regions of space are the modes which 
are occupied by a given number of particles (excitations). 
The entanglement therefore takes the form of non-local 
particle number correlations between the regions and not 
between different Bosons in the gas, or their internal de- 
grees of freedom. This is because even non-interacting 
gases in first quantisation description do have coupling 
between different regions in the second quantisation for- 
mulation and this is what in fact leads to the spatial 
entanglement. A toy model is given in Fig. [1] that helps 
make clear the idea of spatial entanglement as it will be 
central to the rest of this paper. 

In this paper we use two localised probes coupled for a 
short time to two distant regions of a Bose gas to measure 
the spatial entanglement. To measure the entanglement 
correctly we must ensure that the probes act with only 
local operations on the spatial modes. The probes can- 
not couple to the Bosons own degrees of freedom, like 
momentum, as their eigenfunctions do not form a local 
basis in the artificially-defined regions but are meaningful 
only with respect to the entire well. However, the probes 
can couple locally to the number of Bosons in each re- 



FIG. 1: A simple example illustrating why spatial entangle- 
ment must be considered in second quantisation. Both di- 
agrams show a particle in the lowest energy level of an in- 
finite square well. The well has been split into two halves, 
not physically, but conceptually as this allows one to speak 
about regions A and B. Diagram I illustrates the situation 
in first quantisation. In this representation the particle oc- 
cupies the wavefunction as shown and if one wants to speak 
about the state of the particle w.r.t. regions A and B, one 
can say that the particle is in a superposition of A and B. 
However, one cannot speak about entanglement of A and B 
as the Hilbert space of the whole system is a direct sum of 
the Hilbert spaces for A and B and there are no observables 
that we can be measured in regions A or B alone. Diagram 
2 shows the same situation but described in the language of 
second quantisation. In this case the regions A and B become 
the spatial modes and the particle now occupies both the re- 
gions non-locally where the state of the two regions would be 
\iPab) oc |1)a ® |0)s + |0)a ® |l)s- In this representation the 
Hilbert space is a tenosr product of Hilbert spaces of regions 
A and B and we can now measure the number of particles in 
each region. The particle that was represented by a matter 
wave in diagram 1 is still found across the entire well, but the 
superposition of left and right in first quantisation becomes 
particle number entanglement between the regions in second 
quantisation. 



FIG. 2: A confining box of volume V = L contains a Bose 
gas (small spheres) at temperature T. Two two-level systems 
A and B (large darker spheres) interact with two localised 
regions of the gas, Q,a and ils separated by Lab = \r — f '\. 



in diagonal form 



(2) 



one can connect the ODLRO with the eigenvalues 
- the number of particles occupying the single-particle 
states 0£ of a system. For instance a BEC occurs when 
the lowest of the single particle states <f>o is occupied in 
macroscopic way, i.e. Nq/N is finite, where N is the total 
number of particles. We shall see that it is the enormous 
occupation of the ground state mode that accounts for 
the spatial entanglement below the critical temperature 
for condensation (Tc). 



gion, even though the regions are populated non-locally, 
as in this case the number of Bosons in the spatial modes 
form complete bases for the modes. To summarise: be- 
tween the probes we assume no direct interaction and 
any non-local coupling between them must originate in 
the non-local nature of the gas they are coupled to. 



II. THE RELATIONSHIP BETWEEN SPATIAL 
ENTANGLEMENT AND ODLRO 

Let us first discuss the system in question. We take 
a gas of non-interacting Bosons in thermal equilibrium 
described by the temperature density operator 



PG = 



= L e -P(H G -IJ.N) 



z 



(3) 



B. Off-diagonal long-range order (ODLRO) in a 
Bose gas 

In a BEC ODLRO [H, [lj| is present in its simplest 
form [l7| when the one-body reduced density matrix 



Pi 



(1) 



is finite as \f — r'\ — > oo, where ^(r) and \&(r) are 
the creation and annihilation operators of particles at 
point r. The existence of ODLRO testifies that two dis- 
tance points of a BEC have become strongly correlated 
and a new thermodynamic phase emerges, the condensed 
phase. As the one-body density matrix can be expressed 



of the grand canonical ensemble, where Hq is the Hamil- 
tonian of the gas, (3 = 1/ksT where ks is Boltz- 
mann's constant and T is the temperature of the sys- 
tem, Z is the grand partition function and N being 
the number operator for all Bosons. The chemical po- 
tential /i that accounts for the particle number fluctua- 
tions is fixed for any temperature by the implicit relation 

N = E£ i/{f { - E *->* - 1) • 

The gas is contained in a quantisation box of vol- 
ume V = L 3 (Fig. [2]) with periodic boundary conditions 
and the Bosons occupy energy eigenmodes, weighted 
by the temperature, characterised by the plane waves 
4>j:(r) = V~z exp(i/c • r), where fcj = 2irli/Li, i = x,y,z 
and h are the wave quantum numbers in the three spatial 



3 



dimensions. The free Hamiltonian of the non- interacting 
gas is H G = E %a\a& where S g = f v dr cj)^{r)^ (f) 
and are the creation and annihilation operators of 

Bosons in the fc-th mode with energies Er = Ar^-- 



A. Interaction of gas with localised probes 

As we are interested in entanglement between distant 
parts of the gas we need to define two separate regions 
(A and B). This can be achieved in a very natural way 
using the effective volumes of two probe systems VLa and 
ils, which are approximated by spheres of radius R, see 
(Fig.©. 

Here, the probes are two general two-level systems that 
come into contact with the regions A and B for some 
short interaction time tint ■ Probe Sa interacts only with 
region A and probe Sb interacts only with region B. 
The probes start in the unexcited state pp — |00)(00| 
(where |00) = |0)a|0)b) and are therefore initially sepa- 
rable w.r.t. one another. Therefore, as their interaction 
with the gas is local, any entanglement that is found be- 
tween the probes afterwards must have come from the 
particle number entanglement in the g local opera- 
tions alone cannot create entanglement. The interaction 
Hamiltonian 



Hi = Qa&s a + Qbo 



'S B 



(4) 



says that if there are some quanta inside a region (de- 
tected by the operator Qa for region A) , the probe inter- 
acting with that region becomes excited (i.e. ct~s a \0)a = 
|1)a)- A realistic way to implement this interaction is 
discussed in section iHDl 

The full Hamiltonian for the time-evolution consists 
of the free evolution of the gas, Hq and the interac- 
tion between the gas and the probes weighted with some 
coupling r(t) which may be time-dependent, H(t) = 
Hq + T^Hl 

The initial state of the probes and the gas is p(0) = 
Pg ® Pp and if the interaction time ti nt = St is short 
enough one can express the unitary time evolution as 



U = exp ft J dt H(t)j « U(St) « [1 + i(StH G + THi)] 

(5) 

where T — J Q dtT(t). After the interaction the probes 
and the gas are described by the un-normaliscd state 
p{St) = U (St) p(0) W (5t) , which can be reduced to the 
state of the probes, pab, by tracing out the gas degrees 
of freedom. If we consider the limit of very small interac- 
tion time St << I but finite, T << 1 we find the probes 
in the state 



Pab oc tr G [(l +tT(Q A (Tg A +Qb<Ts b )) PG 

®|00)(00| (1 - iT(Q A 4 A + Q B *+ a ))] 
cx |00)(00| + T 2 tr G [(Q A |10) + Q B |01)) . . 



(6) 



x p G (Qa (10| + Q B (01|)] 

-NT tr G [(Qa |10)(00| + Q B |01)<00|)p G ] + h.c. 

After rearranging and normalising we obtain the final 
state of the probes after the interaction written in the 
basis 1 00), |10), |01), |11), 



Pab 



1 

TV 



1 

-iT(Qa) 
-iT(Q B ) 




iT(Q A ) 

r 2 (Qi) 

T 2 (Q a Qb) 




iT{Q B ) 
T 2 (Q a Qb) 

r 2 (g|) 
o 



where (O) denote the expectation values over the gas, 
tr[0 Pg}- Here the normalisation is N = 1 + T 2 ((Qa) + 
(Qb))- Note, that for short interaction times the proba- 
bility of the probes both receiving a kick, represented by 
the right-hand bottom corner of the matrix , the |11) (11 1 
element, representing pab, scales as T 4 which is much 
smaller than the other terms and will be neglected. This 
is an approximation for short interaction times, which we 
shall speak about again later. 
The position observables are 



Qa(b)= df^(f)f Ai B)(r)^(r), 



(8) 



where fA(B)(r) is a real, positive function with support 
only in A(B) that defines the shape of our regions. Here 
we will take /a(b)(' ? ) = C to be equal to a constant inside 
the region and zero outside. The coefficients (Q A ) etc. 
are evaluated in the number basis and arc 



(Qa(i 



drfA(B)(r)pi(r,r), 



(Qa(B)) = J d rf A[B) (r)p l (r 1 r) + 

n 2 J dfdr" f A ( B ) (f) /a(b) {?') g(r, f'\ r, f'), 
(QaQb) = n 2 I 'dfdr' f A (r)f B (r')9(?,?'-A?'). (9) 



where n is the density of Bosons. 

The state of the probes has been related to the reduced 
one-body density operator p\(r,f') = (r)^f(r')) and 
the pair distribution function 

g(r,f';r,r') = n~ 2 {¥(r)^(r')^(r)^(r')) (10) 

of the gas. As the gas is translationally invariant one 
can express p\ and g as functions of only the relative 
distances, r — r' which leads to 



Pi(r-f') = l£(4^)e^'), 

k 



(11) 



g(f-r') = 



1 



Etetata- r a-, r) 



ik-(r-r' 



p,q,k 



where we have used the expresssion ^(r) = J^k ^fc(^)^g- 
Additionally for a non-interacting Bose gas one can then 
write the pair distribution function as the square of the 
one-body density matrix g(r — r') = 1 + (pi^— r')/n) 2 . 
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B. Entanglement of the probes 

Let us take a closer look at the final state of the probes 
Pab ®- The presence of the 1 00} (00 1 accounts for the 
'non-events' where neither of the probes have interacted 
with the gas. This massive 'background' has to be sub- 
tracted to reveal any effects of an interaction and possi- 
bly entanglement arising from it. If we include this term 
in the analysis of entanglement of course we find that 
the probes are in a separable state w.r.t. one another. 
However if we look at the sub-ensemble where the probes 
interacted with the system, we find entanglement. Math- 
ematically we project pab into the subspace spanned by 
the projector P = 1— |00)(00| to the new density matrix 
p' AB . The two probes are now in an entangled state, but 
because we have used a global operation, i.e. non-local 
operation, we can therefore not simply claim that we have 
extracted any entanglement from the gas, we could have 
produced it by our projection. Nevertheless, as the same 
projection has been applied to the probes regardless of 
the temperature of the gas, any temperature dependence 
of the entanglement of the probes comes from a change 
of the spatial entanglement with temperature inside the 
gas. The maximum amount of 'false' entanglement that 
has been generated by this projection will be calculated 
later and will be taken as our zero (background) level. 

We can now investigate the increase in spatial entan- 
glement of the gas from the false background level by 
quantifying the entanglement between probes. By us- 
ing the strategy of inserting two-level probes into the gas 
that pick up the entanglement we have managed to re- 
duce a many-body entanglement problem into a typical 
qub it p roblem and hence we can apply the negativity 
US, LL9|; a standard entanglement measure, to quantify 
the entanglement. The negativity is defined as the sum of 
the negative eigenvalues of a partially transposed density 
matrix [20j | . 

The partial transposed matrix of p'ab-i 



Pab 



(QaQb 

o (Q\) o 

(Q%) 

(QaQb) 

(13) 

reveals one negative eigenvalue and the negativity A/", 
that measures the amount of entanglement in p' AB is then 



N = 



(QaQi 



(Q 2 a) + (Q%) 



(14) 



&+Jn A ^J QB dr' Pl (r-r'f 
2nQ + n 2 Q 2 + J n dr J Q df pi{f- r') 2 

where we have taken the volumes of the two regions to 
be equal, Qa = Qb = & and where /a(b) are taken as 
top-hat functions within the respective volumes and zero 
outside. 

It is clear to see that for a fixed density n and probe 
volume f2, the negativity is affected by the behaviour of 



the reduced one-body density matrix pi(f — r') alone. 
We note again that the presence of ODLRO in a BEC 
depends entirely on the behaviour of the reduced one- 
body density matrix pi(r — r') at large separations. Thus 
the amount of spatial entanglement is directly related to 
the existence of ODLRO. To find out how precisely the 
nature of this dependence is we proceed to calculate the 
negativity explicitly. 



C. Calculation of the negativity 

The presence of ODLRO in a Bose gas is indicative of 
the new condensed phase so our task now is to compute 
the integrals in the negativity N and confirm that there 
is an increase in the amount of entanglement below Tc as 
(Q3J) suggests. There are two kinds of double integrals in 
Af, namely, a cross-over term for the interaction between 
regions A and B 



T AB = \ df dr' Pl (r-r') 2 (15) 



and a on-site term, for instance in the volume Qa 



t aa 
x i 



dr 



dr ' p\(f 



r'f. 



(16) 



For evaluation we use the expression given in pip as 
our starting point and remember that we must treat the 
ground state separately when temperatures below Tq are 
considered. We approximate the sum over momenta in 
equation (|11[) as an integral and evaluate it for low k as 
for temperatures around Tc high momenta modes have 
a infinitesimal occupation. We can then write the one- 
body density operator as 



Pi(f- r') 



exp ( — yjAir{\ 



A- 



( + n 



(17) 

where uq is the ground state density and should only 
be included when the gas is below Tc [21j . The fu- 
gacity is denoted by z = exp(pfkBT) and A is the 
thermal wavelength given by A = \J h 2 /27rmfcsT = 
(2.612 /n) 3 (Tc /T)i , where in the last equality we have 
used an expression for Tc [H|. 

We now evaluate the two integrals above Tc where the 
ground state occupation number no can be neglected. For 
spherical volumes of the probes, Q = ^-R 3 , we obtain 
after some calculation T AA to be 



1 AA (T > Tc) 



2ttV 



'A 4 



2n 2 R 2 



-kR 



(l + 2Ki?)exp(-2 K i?) , (18) 



where k = 2^/47r(l — z)/\. For T AB , the two regions A 
and B are far apart and their radius R shall be much 
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smaller than the distance between them 
tance |r 
in fls 



Thus the dis- 

r ' | in the integral I AB with r in Qa and f 1 
or vice versa is large and can be set to a con- 



stant Lab ~ as it does not change much over the 

comparatively small volume of the individual regions O. 
As we are investigating the relationship of spatial entan- 
glement to ODLRO, the distance Lab between the two 
regions is taken to infinity Lab —> oo, which means that 
the reduced one-body density matrix in the cross-over in- 
tegral (TT5)) vanishes which means the integral is also zero, 
I AB (T > T a ) = 0. 

Below Tc wc must include the ground state density no 
becomes macroscopic and must be included in pi(r — f') 
(fT71 This results in a non-zero cross-over term, X AB (T < 



Tc) 



}f2 , where again the distance Lab — * oo. The 



on-site terms is comprised of three contributions that 
arise from the inclusion of the ground state density in 
the one-body reduced density operator pi(r — r 1 ). The 
full integral 1 AA (T < T c ) was found to be 



I AA (T < T c 



2r>2 



= n^Q, 



l AA (T>T c ) + n l', (19) 



where 



1' 



128tt 2 z 

K 5 A 2 ' 



A-n 2 R z + -n A R i -(A 



-4nR + K 2 R 2 )e- KR }. 

(20) 

Here the ground state density is n = n[\ — (T/Tc)*]. 
We set the ground state energy to zero and use p = 
— /csTln(l + 1/Nq) [2l[ to determine the small correction 
to the chemical potential. 



D. Results and Discussion 

With the integrals computed the focus can shift back 
to the negativity (fT4|) . which should be evaluated for a 
range of temperatures. However, care must be taken to 
chose the parameters carefully as the short interaction 
time approximation, which enabled us to set the coef- 
ficient of the |11)(11| to zero, no longer holds for some 
values of TnQ. In pab 0, the terms that contain the 
highest powers of TnQ are T 2 (Q\) cx T 2 n 2 Q 2 that cor- 
responds to the states |01)(01| and |10)(10| and the ma- 
trix element corresponding to |11)(11| has terms with 
the powers T 4 (QaQbQaQb) oc r 4 n 4 fi 4 . Thus for our 
approximation to be valid the quantity TnQ << 1 must 
stay small so that it is still safe to neglect the |11)(11| 
term from the final state of the probes. In contrast if 
TnQ > 1 even though T « 1, the state |H)(H| would 
participate. 

We would now like to see how the negativity alters 
in the neighbourhood of Tc- To correctly plot the ex- 
pression in (fl4|) we need to weight it with the prob- 
ability that the probes have interacted with the gas 
tv p [PpABP^]- Only if the interaction happened will the 
probes be in the state p' AB with negativity M . The prob- 
ability ti p [P pab P^] for the probes to have interacted 
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Temperature (fraction of T ) 



Density (cm 3) 



FIG. 3: The surface depicts the amount of weighted entangle- 
ment £ (O for T = 2.4 x 10~ 5 , R = 10~ 4 cm and N = 10 6 as 
a function of the temperature of the gas and the particle den- 
sity. Above Tc (Tc — 1 on our scale), the constant amount 
of entanglement between the probes is exactly the maximum 
amount of false entanglement £f (see text). As the temper- 
ature drops below Tc, the amount of entanglement between 
the probes increases significantly above the background level 
£f- This shows that spatial entanglement is directly related 
to the onset of ODLRO below T c . 



with the gas, with (Q 2 A ) = (Q%), is 



tr P [P pab P f } = 



tr p [(l- |00>(00|)pab(1- |00) 

r 2 



N 



tr p [(Oi)|10)(10| + (Q|)|01)(01| 



+«QaQb)|10)<01| + 1i.c.)] 

2r 2 (Q 2 4 ) 



i 



2T 2 (Q 2 A 



(21) 



It is important that the negativity is weighted in this way 
as there is no point talking about entanglement between 
the probes if they have never interacted with the gas. 
The probability is also temperature dependent and the 
resulting quantity is the averaged or weighted entangle- 
ment £ given by 



tr p [P pab P^]M 

T 2 {n 2 Q 2 +l AB ) 
l + 2T 2 {nn + n 2 n 2 +I AA ) 



(22) 



The weighted entanglement is the average entanglement 
between an ensemble of probes where only a certain frac- 
tion of the probes, given by the probability tr p [P pab P }, 
have interacted with the gas. If one could know for cer- 
tain that two probes interacted with the gas then the 
entanglement between them would be given by the nega- 
tivity (|14p . but such entanglement measures are defined 
on ensembles. The plot of £ is shown in Fig. [3] 
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The actual amount of entanglement picked up by the 
probes is extremely small compared to a maximally en- 
tangled state which would have a negativity of one half. 
This is because we have calculated the weighted entan- 
glement and there is a low probability for the probes to 
interact with the gas in the short interaction time. If 
the probes always interacted with the gas, so that the 
contribution in pab of the 1 00) (00 1 term were zero, 
we would end up with a maximally entangled state if 
nil >> 1 at very low temperatures, T — > 0. We can see 
this by looking at the negativity (TH| and noting that 
the integrals X^ A and X' go to zero in these limits. The 
terms with n 2 f2 2 would then dominate and the negativity 
would tend to one half, M —* 1/2. However this is not the 
case in this paper as Tnfl << 1 because we have chosen 
a very short interaction time parametrised by T and very 
small probe volumes f2 « 10 _12 cm 3 . Although note that 
we have a typical condensate density of n = 10 14 cm~ 3 . 

Although the entanglement between the probes is 
small, the amount of weighted entanglement £ behaves 
as expected. It is constant above the critical tempera- 
ture and only sharply increases below the condensation 
temperature Tq where ODLRO is present (see Fig. [3]). 
So indeed, the entanglement between two distant regions 
A and B in the Bose gas, which was picked up by the 
probes, is directly related to the existence of ODLRO in 
the gas, which emerges below Tc- The constant amount 
of entanglement above Tc is artificial and results from the 
global - non-local - operation we made that post-selected 
only those scattering events when the gas and probes in- 
teracted and got rid of the background of non-events. 

We should therefore calculate the amount of false en- 
tanglement £p between the probes due to the projection 
and compare this with the constant amount of entangle- 
ment above Tc- Consider a general, pure and symmetric 
product of the two probes state yS>) = \i/j)a ® W)b where 
\tp) is parametrised by e, = \/l — e|0) + \/e\l). If 
the probes were in this state then they would be sepa- 
rable and the negativity would vanish. However we can 
create entanglement by acting on \^f) globally with the 
projector P = 1 — |00)(00|. To represent the short in- 
teraction time approximation, e << 1 is taken and the 
term where the two probes are both excited, |11), can be 
set to zero. The amount of false entanglement generated 
is £p = e 2 (l — e) 2 /(l — e 2 ), where here the negativity 
has been weighted with the probability of finding the 
probes in the state spanned by the projector P. If we 
choose e = 0.01 so that it is consistent with Tn£l « 0.01 
used in our calculation, the amount of false entangle- 
ment generated by the non-local projection on \^>} is 
£p w 10 -4 . The amount of false entanglement calculated 
here is identical to the constant amount of entanglement 
above Tc for the real state p' AB shown in Fig. [3J namely 
£ = 1.09 ■ 10~ 4 for n = 10 14 cm" 3 . 

This demonstrates that above Tc the entanglement 
between the probes occurs solely because of the global 
projection, but below Tc the increase in entanglement 
between the probes exceeds £f and must come from 



the entanglement of the gas. Interestingly we could use 
£f as the reference level of entanglement, set it to zero 
and use a suitably scaled £ as an order parameter for 
a BEC, which must be zero above Tc and finite below Tc- 

We will now discuss the relationship between our 
work and another paper [22| that also investigates the 
interaction of two probe-type particles with a Bosonic 
field. Firstly we note that the first term of the one 
body reduced density matrix (TTT]) is a Yukawa-type 
interaction potential and the second term is the ground 
state density. We have seen that as the distance between 
the probes goes to infinity the Yukawa interaction term 
goes to zero and only the ground state density is left. 
It is the filling of the ground state that is responsible 
for the ODLRO that is also responsible for the spatial 
entanglement. 

In fact two massive particles interacting via a Bosonic 
field is discussed in [22j and the authors found that below 
the critical temperature for condensation there is an infi- 
nite ranged interaction between the particles, which was 
attributed to the collective occupation of all the Bosons 
in the ground state. We have seen that it is the macro- 
scopic occupation of theground state that is responsible 
for the spatial entanglement in the gas and causes two 
separable probes that interact locally with the gas to be- 
come entangled. By reformulating the problem in terms 
of second quantised spatial modes we have shown that 
two non-local properties of a BEC, spatial entanglement 
and infinite-ranged forces arc equivalent. 

Let us finish this section with a brief discussion of how 
one may implement the probe-gas interaction experimen- 
tally. We need to think of a physical scattering process 
between probe and the gas where a number of bosons 
in a region interacts with an internal degree of freedom 
of the probe. This type of interaction could be realised 
by inserting two probe atoms into the gas that are much 
heavier than the surrounding Bosons. During the short 
interaction time the probes would not move about much 
but an internal degree of freedom, for instance the nu- 
cleus spin, may become excited when hit by a Boson. 

III. CONCLUSION 

In this paper we have presented an operational method 
of investigating entanglement between spatial modes in a 
BEC. Spatial entanglement exists between two distant re- 
gions below the critical temperature for condensation due 
to the massive increase in population of the ground state 
energy level. This is evident in the negativity Eq. (| 14[) . 
due to its dependence on the reduced one-body density 
matrix. As ODLRO is also defined through the behaviour 
of the one-body reduced density matrix and exists only 
below Tc, from the results in this paper we can conclude 
that ODLRO is responsible for the spatial entanglement 
of the gas. 
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